A six-body calculation of the alpha-deuteron radiative capture cross section 
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We have computed the cross section for the process d{a, 7)^Li at the low energies relevant for 
primordial nucleosynthesis and comparison with laboratory data. The final state is a six-body wave 
function generated by the variational Monte Carlo method from the Argonne v\g, and Urbana IX 
potentials, including improved treatment of large-particle-separation behavior. The initial state is 
built up from the a-particle and deuteron ground-state solutions for these potentials, with phe- 
nomenological descriptions of scattering and cluster distortions. The dominant E2 cross section 
is in reasonable agreement with the laboratory data. Including center-of-energy and other small 
corrections, we obtain an El contribution which is larger than the measured contribution at 2 MeV 
by a factor of seven. We calculate explicitly the impulse-approximation Ml contribution, which is 
expected to be very small, and obtain a result consistent with zero. We find little reason to suspect 
that the cross section is large enough to produce significant ''Li in the big bang. 
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I. INTRODUCTION 



Radiative capture of deuterons on alpha particles is the only process by which ^Li is produced in standard primordial 
nucleosynthesis models Because the low-energy cross section for this process is so small (< 10^^ nb), it has long 
been held that no measurable amount of ^Li can be made in big bang nucleosynthesis (BBN) without recourse to 
exotic physics (baryon-inhomogeneous scenarios, hadronically-evaporating black holes, etc.). However, there has been 
new interest in ^Li as a cosmological probe in recent years, for two reasons. First, ^Li is more sensitive to destruction 
in stars than is ^Li. There has consequently been an attempt to set limits on ^Li depletion in halo stars by determining 
how much ^Li they have destroyed over their lifetimes. Second, the sensitivity of searches for ^Li has been increasing 
as a result, so that there are now two claimed detections in metal-poor halo stars, and two more detections in disk 
stars [||-^. The explanation of these data in terms of chemical-evolution models remains an open question [^-^. 
They are presently at a level exceeding even optimistic estimates of how much ^Li could have been made in standard 
BBN and the observed ^Li is believed to have been created in interactions between cosmic rays and the interstellar 
medium. However, it remains interesting from the point of view of understanding these observations to remove as 
many remaining uncertainties as possible from the standard scenario. The main such uncertainty arising in the 
standard BBN calculation is the cross section for d{a^ 7)^Li- 

At the same time, electroweak processes in light nuclei are interesting from the point of view of few-body nuclear 
physics. The advent of realistic two- and three-nucleon potentials has enabled calculation of nuclear wave functions and 
energy observables of systems with up to A = 8 ||] . Calculation of electroweak observables allows both comparison 
of these wave functions with experiment {e.g., electron scattering form factors and estimates of effects not 

observable in the laboratory [e.g., weak capture cross sections The radiative capture of deuterons on alpha 

particles provides both kinds of opportunities: there are data on the direct process at 700 keV (all energies are in 
center of mass) and above 14 1, and there are indirect data from Coulomb breakup experiments corresponding to 
70-400 keV There is also a need for extrapolation to lower energies for application to nucleosynthesis. 

The ad capture problem is made even more interesting by the fact that S- and P-wave captures are strongly 
inhibited by quasi-orthogonality between the initial and final states and by an isospin selection rule, respectively. 
As a result, the dominant process in all experiments performed to date has been electric quadrupole (i?2) capture 
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1. Sensitivity function of Ref. showing the change in ^Li BBN yields resulting from small changes in the d{a, 7)^Li 




cross section over narrow energy bins at the specified energies. 



from Z?-wave scattering states. The small remaining El contribution from P-wave initial states has been observed 
at about 2 MeV, but its magnitude has not been successfully explained by theoretical treatments; it is generally 
expected to contribute half of the cross section at 100 keV. The 5-wave capture induced by Ml has been neglected 
in most calculations because of the quasi-orthogonality mentioned above, which makes the associated matrix element 
identically zero in two-body treatments of the process. The energy dependences of the various capture mechanisms 
{E2, El, Ml) are such that even El and Ml captures with small amplitudes may become important at low (< 200 
keV) energies. Low-energy behavior is particularly important for standard BBN: the primordial ^Li yield is only 
sensitive to the capture cross section between 20 and 200 keV, with the strongest sensitivity at 60 keV. This is 
demonstrated in Fig. |l], where we show the fractional change in ^Li yields resulting from changing the cross section 
over narrow bins in energy. Since it is the energy integral of this quantity which determines the sensitivity, and it is 
shown on a logarithmic axis, the sensitivity function ge is defined to include a factor of energy so that relative areas 
may be accurately gauged. (A more extensive discussion of these functions appears in Ref. The ^Li yield is 

directly proportional to the c?(a, 7)^Li cross section at the energies indicated by this function, so that an increase in 
the cross section by a given factor over the whole energy range produces an increase of the ^Li yield by the same 
factor. 

We have carried out a calculation of the alpha-deuteron radiative capture cross section, treating it as a six-body 
problem. The remainder of this paper describes this calculation, and it is organized as follows: In Section ||, we 
describe wave functions used to compute the £i(a,7)^Li cross section. In Section |l|, we describe the electromagnetic 
current operators and the methods used to calculate their matrix elements. In Section |l^, we describe our results for 
the cross section, thermal reaction rates and nucleosynthesis yields, and in Section ^ we summarize our results and 
discuss their implications. 



II. WAVE FUNCTIONS 



A. ^H, *He, and "^Li ground states 



The wave functions IV'™'*)' IV'a): ^iid ipiT) used in our calculation are the ground states derived from the Argonne 
Vis two-nucleon and Urbana IX |l8] three-nucleon potentials, henceforth denoted as the AV18/UIX model. The 
deuteron wave function is a direct numerical solution, while the variational Monte Carlo technique described in 
Refs. [^-21 is used to generate the *He and ^Li wave functions. Here and mg denote spin orientation. 

The variational trial functions for light nuclei are constructed from correlation operators acting on a Jastrow wave 
function: 



i<j<k 



TNI 
ijk 
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where Uij and U^^^ are two- and three-body correlation operators that include significant spin and isospin dependence 
and 5 is a symmetrization operator, needed because the Uij do not commute. For ^He, the Jastrow part takes a 
relatively simple form: 



i*j> = n f^i" n /(^»j)i*o(oooo)i234) 

i<j<k<4 i<j<4: 



(2) 



where f{rij) and fijk are pair and triplet functions of relative position only, and (f>Q(0000) is a determinant in the 
spin-isospin space of the four particles. 

For ^Li the Jastrow part has a considerably more complicated structure due to the need to place the fifth and sixth 
particles in the p-shell: 



n f^IkUf^ n Mr^J)llfsp{rk5)fspirk6)fpp{r,6) 
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The A is an antisymmetrization operator over all partitions of the six particles into groups of four plus two. The 
central pair and triplet correlations fxyi^ij) and f^^^jf have been generalized, with the xyz denoting whether the 
particles are in the s- or p-shell. The wave function \^q{LS J MTT^)) has orbital angular momentum L and spin S 
coupled to total angular momentum J, projection M, isospin T, and charge state Ts, and is explicitly written as 

|$6(^^^MrT3)i234:56> = (00 00) 12340^^ (iJas)^^'^ (i?a6) 

^. ' J JM 

x[,y,{^tMht3)]Tn) ■ (4) 



Particles 1-4 are placed in the s-shell core with only spin-isospin degrees of freedom, while particles 5-6 are placed in 
p-wave orbitals 4>p^{Rak) that are functions of the distance between the center of mass of the core and the particle. 
Different amplitudes Pls are mixed to obtain an optimal wave function; for the J'^, T = 1"*", ground state of ^Li the 
p-shell can have /3oi, /32i, and (3iq terms. 

The two-body correlation operator Uij is defined as: 



E 

p=2,6 



n f^kiT^^k^Yjk) 



UpinjMj , (5) 

p=2,6 k=jti,j 

where the 0^~'^'^ = Ti-Tj, cTi crj, cri-a-jTi-Tj, Sij, and SijTi-Tj, and the f^^f. is an operator-independent three-body 
correlation. The six radial functions fssir) and Up=2s{r) are obtained from two-body Euler-Lagrange equations with 
variational parameters as discussed in detail in Ref . [|l9[ . Here we take them to be the same as in ^He, except that the 

Up=2s{'>') are forced to go to zero at large distance by multiplying in a cutoff factor, 1 + exp[— i?„/a„] / 1 + exp[(r — 



, with Ru and a„ as variational parameters. The fap correlation is constructed to be similar to fss for small 



separations, but goes to a constant of order unity at large distances: 

h 



fsp{r) 



•'sp 



''sp 



1 + exp[(r - Rsp)/asp] 



fss{r) + Cspil ~ exp[~{r/dsp) ]) 



where asp, hsp, etc. are additional variational parameters. The fpp(r) correlation is given by: 

fpp{r) = Udir)/[1 - 3u2(r) + u^ir) - 3M4(r)]r , 



(6) 



(7) 



where Ud(r) is the S'-wave part of the exact deuteron wave function. Then Uij acting on the pair of nucleons in the 
p-shell will regenerate the correct deuteron S'-wave and an effective D-wave Wd{r) — fpp{r)[u^{r) — 3ug(r)]-\/8?' from 
the tensor components. While the D-wave is not exact, we find that the resulting effective deuteron energy is —2.16 
MeV, quite close to the correct value. 
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These choices for f^s, fsp, fpp, and Up=2, 
core, the overall wave function factorizes as: 



guarantee that when the two p-shell particles are far from the s-shell 



(8) 



where ipa is the variational ^He wave function and tpd is (almost) the exact deuteron wave function. 

The single-particle functions (pp^jMan) describe correlations between the s-shell core and the p-shell nucleons, and 
have been taken in previous work to be solutions of a radial Schrodinger equation for a Woods-Saxon potential 
and unit angular momentum, with energy and Woods-Saxon parameters determined variationally. It is important 
for low-energy radiative captures that these functions reproduce faithfully the large-separation behavior of the wave 
function, because the matrix elements receive large contributions at cluster separations greater than 10 fni. In fact, 
below 400 keV, more than 15% of the electric quadrupole operator comes from cluster separations beyond 30 fm. We 
have therefore modified the ^Li wave function for the capture calculation to enforce cluster-like behavior when the 
two p-shell nucleons are both far from the s-shell core. 

In general, for light p-shell nuclei with an asymptotic two-cluster structure, such as ad in ^Li or at in ''Li, we want 
the large separation behavior to be 



if 



(9) 



where Wkm{'2jr) is the Whittaker function for bound-state wave functions in a Coulomb potential (see below) and n 
is the number of p-shell nucleons. We achieve this by solving the equation 




r</)f (r) = 0, 



(10) 



with i — 1, flu the reduced mass of one nucleon against four, and V{r) a parametrized Woods-Saxon potential plus 
Coulomb term: 



y(r) = 



Vn 



2{Z-2) : 



1 + cxp[(r - i?o)/ao] 



-F{r) 



(11) 



Here Vq, Rq, and uq are variational parameters, {Z — 2)/n is the average charge of a p-shell nucleon, and F(r) is a 
form factor obtained by folding a and proton charge distributions together. The A(r) is a Lagrange multiplier that 
enforces the asymptotic behavior at large r, but is cut off at small r by means of a variational parameter cq: 



A(r) = A(r)[l-exp(-(r/co)2) 



The A(r) is given by 



(12) 
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where ml is directly related to the Whittaker function: 



(14) 



Here 7^ = 2/i4„_B4„/?i , with /i4„ and B^n the appropriate two-cluster effective mass and binding energy, k — 
-2{Z - 2)e2^4ri/?i^7, and m i + i. 

For ^Li, ,642 = 1.47 MeV and i = 0, 2 corresponding to the asymptotic S- and Z?- waves of the ^Li ground state, 
or amplitudes /3oi and /32i in Eq.(H). There is no asymptotic ad cluster corresponding to the L = 1 amplitude (3io. 
Nevertheless, the energy of the ^Li ground state is variationally improved « 0.1 MeV by including such a component 
in the wave function, so we set the asymptotic behavior of 4>p^{r) to a binding energy of 3.70 MeV, which is the 
threshold for ^Li — > ^He + p + n breakup. 

In Refs. the ///^'^ three-body correlation of Eq.(|) was a valuable and inexpensive improvement to the trial 

function, but no f^^j^ or f^^^ correlations could be found that were of any benefit. However, in the present work we 
find the correlation 



/n56 = 1 + <li[fss{r56)/fpp{r56) ~ l]exp[-g2 (^-nS + rns)] 



(15) 
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TABLE I. Calculated VMC, GFMC, and experimental energies and point proton rms radii of ^H, ^He and ^Li. Numbers in 
parentheses are Monte Carlo statistical errors, or experimental uncertainties. 



INucieus 


Ubservable 


VMC 




Cr MC 


Experiment 




E 


-2.2246 






-2.2246 






1.967 






1.953(3) 


*He 


E 


-26.89(3) 


-27.78(3) 


-28.34(4) 


-28.30 






1.48(1) 


1.47(1) 


1.45(1) 


1.48(1) 


E 


-27.19(5) 


-28.21(4) 


-31.15(11) 


-31.99 






2.46(2) 


2.46(2) 


2.57(1) 


2.43(4) 



with qi 2 as variational parameters, to be very useful. It effectively alters the central pair correlation between the 
two p-shell nucleons from their asymptotic, deuteron-like form to be more like the pair correlations within the s-shell 
when the two particles are close to the core. This correlation improves the binding energy by « 0.5 MeV. 

In Rcfs. §,|ll we reported energies for the trial function '^t of Eq.(|l|), and for a more sophisticated variational wave 
function, which adds two-body spin-orbit and three-body spin- and isospin-dependent correlation operators. The 
gives improved binding compared to \E't in both ''He and ^Li, but is significantly more expensive to construct 
because of the numerical derivatives required for the spin-orbit correlations. In the case of an energy calculation the 
derivatives are also needed for the evaluation of L-dependent terms in AV18, so the cost is only a factor of two in 
computation. However for the evaluation of other expectation values the relative cost increase is ~ 6A. Thus in the 
present work we choose to use ^'t for non-energy evaluations; this proved quite adequate in our studies of ^Li form 
factors 0. 

The variational Monte Carlo (VMC) energies and point proton rms radii obtained with and are shown 
in Table | along with the results of essentially exact Green's function Monte Carlo (GFMC) calculations and the 
experimental values. We note that the underbinding of ^Li in the GFMC calculation is the fault of the AV18/UIX 
model and not the many-body method; it can be corrected by the introduction of more sophisticated three-nucleon 
potentials [ p2[ . 

The present variational trial functions with the imposed proper asymptotic behavior in fact give slightly better 
energies than the older shell- model- like correlations of Refs. [ppl|]. Unfortunately, because the variational ^Li energy 
is not below that of separated alpha and deuteron clusters, it is possible to lower the energy significantly by making 
the wave function more and more diffuse. Therefore, the variational parameters were constrained to give a reasonable 
rms radius, as seen in Table ||. Although the variational ^Li wave function remains unbound relative to alpha-deuteron 
breakup, it has the correct energy at large particle separations, as shown in Fig. |^. There we show the local energy, 
i?(R), which is the energy for a given spatial configuration R = (ri, . . . ,rg), binned with its Monte Carlo statistical 
variance as a function of the sum R = |ri| of the particle distances from the center of mass. The average energy 
is shown by the solid line, while the dashed line gives the experimental binding. At the bottom of the figure is 
a histogram of the number of samples to show their distribution with R. We observe that the energy in compact 
configurations is clearly too small, but for i? > 18 fm the energies scatter about the experimental binding. 

The asymptotic two-cluster behavior of our six-body wave function can be studied by computing the two-cluster ad 
distribution function, {Aipa^'T'^ ^''^ad \ described in Ref. ||2^. It can be expressed in terms of Clebsch-Gordan 

factors, YLMi^ad), and two radial functions RhiTad) for the S- and _D-waves, which are plotted in Fig. |. The Rl can 
be used to extract the asymptotic normalization constants Co and C2, by taking their ratio with the corresponding 
Whittaker functions, and the asymptotic D/S ratio "q. We find that Co does not become asymptotic until Tad ^ 8 fm, 
where we obtain Co — 2.28 ± 0.02. However, r] is already asymptotic for r^d ^ 4 fm, with the value —0.026 ± 0.001. 

This value of Cq is at the lower end of the range of values (2.3-2.9) used in a number of other studies [p4[-p6[. 
However, the value of 77 is consistent with the recent analysis of ^Li -|- ''He elastic scattering which obtains 
—0.025 ± 0.006 ± 0.010. This is a big improvement over our earlier value of —0.07 ± 0.02 reported in Ref. |23], which 
was obtained with the shell-model-like wave function. However, our quadrupole moment Q — — 0.71±0.08 fiii^ remains 
an order of magnitude too large compared to the experimental value of —0.08 fm^. The integrated D'^'^ parameter 
value is -0.27 ±0.01 fm^. 

The spectroscopic factor is 0.86, in good agreement with the Robertson value of 0.85 ±0.04 The spectroscopic 
factor is not very sensitive to the asymptotic behavior of the wave function, and differs little from our earlier value of 
0.84 obtained with the shell-model-like wave function. We note that our shell-model-like wave functions for ^Li and 
^He provided an excellent prediction for the spectroscopic factors observed in ^Li(e,e'p) scattering [ p7[ . 
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FIG. 2. Monte Carlo local energy of the ^Li variational ground state binned as a function of the sum of the magnitudes of 
the particle coordinates, R. The horizontal solid line is the average variational energy and the horizontal dashed line indicates 
the experimental binding. The curve superimposed on the bottom is a histogram of the location of the 50 000 samples of ^'y 
used in the calculation. 
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B. ad initial state 



The initial alpha-deuteron scattering wave function \ipad] LSJM) having orbital angular momentum L and channel 
spin S (=1) coupled to total JM, is expressed as 



where curly braces indicate angular momentum coupling, A antisymmetrizes between clusters, and Vad denotes the 
separation between the centers of mass of the a and d clusters. 

Our scattering wave function is built from the alpha and deuteron ground-state wave fuctions, V'a and -0™'* (where 
nid denotes the deuteron spin projection), from a nucleon pair correlation operator Gy , and from an intercluster scalar 
correlation 4'adij'ad)- Gij is the identity operator if particles i and j belong to the same cluster (alpha or deuteron). 
Otherwise, it is a pair correlation operator derived to describe nuclear matter |^^ , and it contains both central and 
non-central terms. It takes into account, in an approximate way, the distortions introduced in each cluster by the 
interactions with the particles in the other cluster, and becomes the identity operator as particle separations become 
large. 

The intercluster correlation 4'adi'>'ad) is a two-body radial wave function describing alpha-deuteron scattering as the 
scattering of two point particles. It is generated from a Schrodinger equation, with the potential chosen to reproduce 
the results of phase-shift analyses of alpha-deuteron scattering data in the two-body model. Several such potentials 
exist in the literature, and we saw no need to re-fit the data. The potentials which we have applied to this study 
are the potentials of Kukulin and Pomerantsev, Refs. [^,0 (hereafter, "KP I" and "KP II", respectively), and the 
potential of Langanke, Ref. Because the potential of Langanke was produced for the slightly different problem 
of describing the ad system in the resonating group method, it required modification to match our purposes. This 
consisted of a small adjustment to the Woods-Saxon radius to reproduce the energy of the S"*" resonance derived from 
electron scattering, 711 keV. The most important difference between these potentials, which give very similar cross 
sections in our calculations, is that the KP potentials were fitted to the phase shifts in the odd and even partial waves 
separately, to produce the expected parity-dependent potential. The work by Langanke was concerned only with even 
partial waves, and did not fit the odd partial waves, although it provides a fair qualitative description of the P-wave 
phase shifts at low energy. 

The potentials are very similar to each other in form. All of the potentials used contain both central and spin-orbit 
terms, and the spin-orbit force in the even partial waves is well-constrained by the spacing of the D-wave scattering 
resonances. While a tensor component could in principle be derived to reproduce the ground-state quadrupole moment 
(assuming total alpha-deuteron cluster parentage of the ground state), its effect on the scattering data is too small to 
support deriving it from the scattering phase shifts. (A more ambitious approach to derive a tensor interaction would 
involve computing energies at fixed separation between the a and d clusters and varying the deuteron orientation; 
however, we have not pursued such an approach here.) Coupling to other {e.g., cluster breakup) channels has also 
been neglected; this is well-justified below the ^He + p and ^Li -|- n thresholds at 3.12 and 4.20 MeV, respectively. 
We note that although the KP potentials and the Langanke potential all reproduce the scattering phase shifts very 
well, and they produce very similar ground-state wave functions at small cluster separation, the ground-state binding 
of the KP II potential is too large by almost a factor of two. All of the potentials we used to generate (pad are 
deep potentials with a forbidden zero-node state. They produce the expected one-node structure of the ground-state 
alpha-deuteron wave function required by the Pauli principle, as illustrated by our Rl functions in FigjH and enforce 
(by orthogonality to the forbidden state) the most important consequences of the Pauli principle for (pad- 

A complication arises from the phenomenological treatment of the initial state. In reality, the ground-state and 
5-wave scattering-state wave functions are orthogonal, because they are different eigenstates of the same ( J'^ = 1+) 
potential. Because our calculation does not generate them explicitly as such, our ground state and S'-wave scattering 
states are not necessarily orthogonal. For purposes of our capture calculation, this has the effect of generating 
spurious contributions to the Ml transition operator. For S'-wave captures, we have modified the initial state to 
ensure orthogonality with the ground state. This is done by subtracting from the initial state, Eq.(^6|), its projection 
onto a complete set of ground states. 



and computing Ml transitions from the adjusted state. This procedure reduced the Ml contribution to our computed 
cross sections by more than a factor of 10, roughly to the level of Monte Carlo sampling noise. Further evidence that 




(16) 




(17) 
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the remaining contribution is spurious is its strong dependence on the phenomenological radial correlation in the 
scattering state, discussed below. 



III. CROSS SECTION AND TRANSITION OPERATORS 



The cross section for ad radiative capture at cm. energy Ec.m. can be written as 



3 -Wi-oi 1 + 9/ 



E 

LSJl L 



LSJ ( 



(18) 



where a is the fine structure constant (a — e^; we take fi 



1 in this Section), is the ad relative velocity, 



and Ef^'^{q) and M^^"^ [q) are the reduced matrix elements (RMEs) of the electric and magnetic multipole operators 
connecting the ad scattering state in channel LSJ to the ^Li ground state having J^,T = 1~^,0. The cm. energy of 
the emitted photon is given by 



mu 



1 + W IH [md + TOq - mu + Ec.m.) 

mu 



md + ma - mu + Ec 



(19) 



where rud, m^, and mu arc the rest masses of deuteron, '^He, and ^Li, respectively. The astrophysical S'-factor is 
then related to the cross section via 



S{Ec.„i.) = Ec.m. <7{Ec.m.) exp{4:n a/vrci) ■ (20) 

The Ml, Af2, El, E2, and Ei transitions involving ad scattering states with relative orbital angular momentum up 
to L = 2 have been retained in the evaluation of the cross section. Initial states with L — i and L — A were examined 
in the early phases of this work, but proved to be very small, and were not retained in the final calculation. 

Since the energies of interest in the present study are below 4 MeV, and consequently qRu < 0.05 [Ru is the ^Li 
point rms radius) , the long- wavelength aproximation (LWA) of the multipole operators would naively be expected to 
be adequate to compute the associated RMEs. While the LWA is indeed sufhcient for E2 transitions, which dominate 
the cross section at energies Ec.m. > 400 keV, it is however inaccurate to compute the RMEs of Ml and, particularly, 
El transitions in the LWA since these are suppressed. (This fact has already been alluded to in the introduction). 
It is therefore necessary to include higher order terms, beyond the leading one which normally is taken to define the 
LWA, in the expansion in powers of q, so-called "retardation terms" . It is useful to review how these corrections arise 
for El transitions. We should point out that some of these same issues have recently been discussed in Ref. [Q, in 
the context of a calculation of El strength in dp radiative capture at energies below 100 keV. 



A. The El operators 



In principle, the El RMEs are obtained from 



= - / <ixj(x) • V X 3,(5a:)Yj;(i) , (22) 

where j(x) is the nuclear current density operator, ji{qx) is the spherical Bessel function of order one, and Y]^j^(x) are 
vector spherical harmonic functions. For ease of presentation, the factor ^/i/ {JM, lAjlwe) occurring in the definition 
of the RMEs will be suppressed in the following discussion. By expanding the Bessel function in powers of q and after 
standard manipulations p^ , the Ei operator can be expressed as 

£;iA(g) ~ -El A (9; LWAl) + E^xiq; LWA2) + E^,^{q- LWA3) , (23) 

where 
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£;iA(g;LWAl) 



V2 
3 



J dxxriA(x)p(x) 



(24) 



Eixiq; LWA2) - ^ y" dx X Fi A (i) X • j (x) 



(25) 



£;iA(g;LWA3) 



V2q' 
15 



H, J dxa;3riA(x) p(x) 



(26) 



Here the continuity equation has been used to relate V ■ j(x) occurring in £^ia((J'; LWAl) and £'iA(g;LWA3) to the 
commutator —i[H, p(x)], where p(x) is the charge density operator. Evaluating the RMEs of these operators leads 
to 



E^^'-'iq) ~ E^^'-'ig; LWAl) + E^'"' {q; LWA2) + E^''-'{q; LWA3) 



(27) 



with 



"^(^ul / dxxriA(i)p(x)|V„,;i5JM) 



£;f^^(q;LWAl) 



(28) 



i?f^-^(<Z;LWA2) = ^(^Lil / dxa;yiA(x)x.j(x)|Vad;L5JM) , 



(29) 



i?f^^(g;LWA3) = y dxa;3riA(i)p(x)|^„rf; L^JAf) , 



(30) 



where terms up to order q'^ have been retained, and it has been assumed that the initial and final VMC wave functions 
are exact eigenfunctions of the Hamiltonian - which, incidentally, is clearly not the case, see the previous Section - 
so that [H , p(x)] — > —qp{x} in the matrix element, ignoring the kinetic energy of the recoiling ^Li. 

The nuclear electromagnetic charge and current density operators have one- and two-body components. 



p(x) = ^pW(x) + ^pif(x) 



i<j 



j(x) = Eji''w+Ej?w 



(31) 
(32) 



The one-body terms pl^'' and have the standard expressions |34| obtained from a non-relativistic reduction of the 
covariant single- nuclcon current. The charge density is written as 



PrHx)=pW (x)+pW (x) 



with 



pInrW = e»'5(x-r,) 



(33) 
(34) 



a.xpi, (5(x - r^) 



(35) 



-I + 



where the Darwin- Foldy relativistic correction has been neglected in the expression for p| since it gives no contri- 



bution to Eq. (28) (see below). Hereafter, [• • • , ■ ■ ■]+ denotes the anticommutator. The one-body current density is 
expressed as 



Pi , (5(x - Tj) 



2m 



Ml CTi X Va;(5(x - Ti 



(36) 
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The following definitions have been introduced: 



(37) 
(38) 



where {ii^— 0.88 n.m.) and /i^ —4.706 n.m.) are the isoscalar and isovector combinations of the proton and 
neutron magnetic moments. 

The two-body charge and current density operators have pion range terms and additional, short-range terms 
due to heavy meson exchanges [Q. However, the pion-exchange current density is an isovector operator, and its 
matrix element vanishes identically in the radiative ad capture, while the heavy-meson exchange (charge and current) 
operators are expected to give negligible contributions. Hence, in the present study we only retain the charge density 
operator associated with pion exchange. It is explicitly given by 



2 m mi 



(Tj ■ VijCTi ■ Vx (5(x - Vi 



(39) 



where TOtt is the charged pion mass, /tt the pion-nucleon coupling constant {f'^/Air — 0.075), 
separation of nucleons i and j, and the function lTi{r) is defined as 



Tj is the 



-A„r 



i ( 1 - ) A.re 



-A„r 



(40) 



The prime denotes differentiation with respect to its argument. For the parameter A^, we have taken the value 1.05 
GeV. 

We now consider, in turn, the various contributions LWAl, LWA2, and LWA3 to the Ei RME. Insertion of the NR, 
RC, and 7r-exchange charge density operators in Eq. (pq) allows us to write correspondingly 



E^^'^iq; LWAl) ~ Ei^^'-'iq; LWAc) + E^'^'' (q; LWAb) + ■^•^ (q; LWAtt) 



where 



^f'^-^(g;LWAc) 



(41) 



(42) 



i?P^(g;LWAb) = ^ ^ p,),]^;^,; LSJM) 



3 V 47r" 



(43) 



E^^^iq; LWA.) = E • I^in,) 



2 m ml 



\ij^d;LSJM) 



where {a x p)x and ai \ denote spherical components. 

Next, for the terms of order q^ and q"^, Eqs. (p9|) and (|30|), we find: 



i;f*''(q;LWA2) 



3\/2 



2m 



Pi-, nYixirt) 
3 



2 m V 47r 



(r^ X (Ti)x\tpad;LSJM) , 



V2t 



E'^s^q- LWA3) = E ri,(f,)|Vod; LSJM) 



(44) 



(45) 
(46) 
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where in the last equation the contributions from relativistic corrections and pion-exchange have been neglected. 

The leading contribution '^'^(q; LWAc) to the Ei RME vanishes in ad capture, since the initial and final states 
have zero isospin, and the corresponding wave functions do not depend on the cm. position: 

E^'^-'iUNAc) « {^'^{>\Y,r^,^\^pc.d;LSJM) = , (47) 

i 

where the coordinates are defined relative to that of the cm.. 

There are also more subtle corrections, which can potentially lead to additional, (relatively) significant Ei contri- 
butions. One such correction is the possibility of isospin admixtures with T > in the ^Li and ad states, originating 
from the isospin-symmetry-breaking interactions present in the AV18/UIX model. These T > components are 
ignored in the VMC wave functions used here. 

An additional correction - which we take into account - arises from a relativistic correction to the center of mass, 
and makes the LWAc contribution nonzero. This correction arises because translation-invariant wave functions require 
not center-of-mass but center-of-energy coordinates. Such coordinates constitute the relativistic analog of the center- 
of-mass coordinates appropriate to nonrelativistic problems. The distinction is usually not crucial at low energies, 
but it matters here because the LWAc term would otherwise be zero. Since the coordinates in which we compute 
wave functions are not center-of-energy coordinates, the dipole moment 



E 



e^r^, (48) 



which enters Eq. (p7|), should be replaced by the expression 

d = ^e,(r,-rcis) (49) 

]<^ircE , (50) 



where 
and 

Ei = m 



rcE = ^ , (51) 



Here, and Vij^ are the two- and three-nucleon potentials. We have thus replaced the particle masses in the 
expression for center of mass by their total energies (rest energy plus kinetic and potential energies). Since the 
nucleon masses contribute most of the nucleons' relativistic energies, the center of mass is only slightly different from 
the center of energy. (The rms radius of the center of energy for our ^Li variational ground state is (9.1 ± 2.7) x 10^'^ 
fm.) We note that this correction has been taken into account in the past (in two-body models) by putting the 
measured alpha and deuteron masses into calculations, rather than 4m and 2m p!3| , p^ , p5| , p5p^ ] . These differ by the 
mass defects of the alpha particle and deuteron. 

The El cross sections arising from the center-of-energy correction are larger than the El cross sections measured 
in Ref. [|l3| by about a factor of six. There is a significant (~ 30%) LWAb spin-orbit contribution which becomes 
more important with increasing energy; the next-largest contributions are from the LWA3 and LWA2 operators. All 
of these terms except LWAb are too small in our calculation to account for the observed transition by more than a 
factor of ten, but LWA3 is responsible for about 20% of the total S factor at zero energy. We have also included the 
contribution labelled LWAtt in the present study, and we find it to be much smaller than LWAb. Note that although 
the LWAb term is large enough to account for the magnitude of the El data alone, it is of the wrong sign to account 
for the asymmetry from which the El strength is extracted. 



B. The Monte Carlo Calculation 



We calculated matrix elements by modifying the variational Monte Carlo code described in Ref. to perform 
Monte Carlo integrations of matrix elements between the scattering and ground states described above. The Monte 



11 



10 20 30 40 10 20 30 40 




I — 1 cluster separation (fm) 

FIG. 4. Upper panel: The integrand of Eq.(M) for E2 transitions from the Di initial state and El transitions from the 
P2 initial state. The upper panels are for Ec.m. = 50 keV, while the lower panels are for Ec.m. = 1 MeV. "Noisy" curves 
are computed solely from the Monte Carlo integration (and their scatter indicates sampling uncertainty), while smooth curves 
show the expected long-range behavior of the operators, with normalization fitted to the Monte Carlo results. 



Carlo algorithm used was the Metropolis algorithm, and we tried several weighting functions, all based on the ^Li 
ground state. The final calculation used the weighting function 



M^(R) = ^«;«(R)|<r(R)) (53) 

(where the bra-ket product here indicates summation over spin and isospin degrees of freedom for a given spatial 
configuration R = (ri,r2, . . . jrg) of the particles; in all other cases, it also indicates integration over space coordi- 
nates). This weighting function was chosen to reflect approximately the expected behavior of the matrix element, 
and to obtain good sampling both at small cluster separation and out to 40 fm cluster separation. In the low-energy 
calculation, 3% of the El matrix element and 15% of the E2 matrix element come from the region beyond 30 fm 
cluster separation. 

Because all of the energy dependence is contained in the relative wave function (j>ad{rad) and the transition operators, 
the matrix element for a given scattering partial wave can be re-written as 

T['^'iq)=J^^dxx'(l)^dixMn'\Tixiq)Als{x~r^d)^^^^^^ \ , (54) 

I '■^ J LSJM 

using a technique from Ref. Here Tex denotes any of the Ei and Me operators. The partial- wave expansion of 
the current operator contains the photon energy only in multiplicative factors on each term of the expansion. The 
integration over all coordinates except x can therefore be calculated once for each partial wave by the Monte Carlo 
code, and the result can then be used to compute the full integral for as many energies as desired by recomputing 
(j)ad only, greatly reducing the amount of computation. 

The calculation sampled 2 000 000 Monte Carlo configurations, summing over all 15 partitions of the six particles 
into clusters at each step. Typical results for the integrand in Eq.(|5^) at low energy are shown for two transitions, 
at two different energies, in Fig. ^ Calculations at energies above about 1 MeV had integrands that were essentially 
zero beyond 30 fm cluster separation, and therefore suffered little from limited sampling. 

We also performed a second set of calculations, this time computing functions which had the expected asymptotic 
forms of the operator densities as a function of x at large cluster separations, and normalizing them to the Monte- 
Carlo sampled operator densities between 12 and 22 fm. These asymptotic forms, given by multiplying Eq. (^ by 
powers of rad appropriate for each operator, were then used instead of the Monte Carlo operator densities at cluster 
separations greater than 12 fm, in order to check that our results were not affected by the effective cutoff in sampling 
beyond 45 fm (due to very small PF(R)) or by poor sampling elsewhere in the tails of the wave function. For all of the 
dominant reaction mechanisms, this gave essentially the same result as the explicit calculation. The smooth curves 
of Fig. ^ were computed in this way. 
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5. The total S-factor calculated from the Langanke (solid curve) and KP II potentials (dashed curve), with the data 
(A), (x), and ^ (o). The indirect data of Ref. || were taken from their graphical presentation in P8|. 



120 
100 
80 

I 

^ 60 
40 
20 






it : 
\ - 











0.68 0.7 0.72 0.74 

E.„,. (MeV) 

FIG. 6. The calculated total cross section at the 711 keV resonance, compared with the data of Ref. |l^. The solid curve 
was computed from the Langanke potential, while the dashed curve was computed from the KP II potential. 



IV. RESULTS 



Our main result is the total capture cross section as a function of energy, shown in comparison with laboratory 
data in Fig. || as an astrophysical S'-factor. Note the good agreement with the capture data below 3 MeV, reflecting a 
reasonable value of the alpha-deuteron asymptotic normalization coefficient of the ^Li ground state, Co, and accurate 
reproduction of the scattering phase shifts. The disagreement with the data above 3 MeV is a generic feature of direct- 
capture models, and is usually attributed to neglected couplings to breakup channels, and to channels with nonzero 
isospin in the scattering state (such as the J'^, T = 0+, 1 state at 2.09 MeV). As shown in Fig. H, the behavior at the 
3+ resonance is nicely reproduced in location, width, and amplitude, simply by constraining the cluster potentials in 
the D waves to reproduce the experimental phase shifts. The failure to match the energy dependence of the cross 
sections inferred below 500 keV from Coulomb-breakup experiments is shared by all other direct capture calculations, 
and could reflect either something that is not understood about the experiments, or something that is missing from 
the models. Our result with more detailed wave functions and operators than previously considered may be taken as 
a further indication that previous models were correct in this region. 

In Figs. 1^ and ^, we show the breakdown of the cross section into contributions from different multipolarities and 
incoming partial waves, again as S factors. Fig. ^ shows results of two calculations, one from the Langanke potential 
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12 3 

E^,,, (MeV) 

FIG. 7. Contributions of the more important reaction multipolarities and incoming partial waves to the S factor. Curves 
marked with the symbols x were computed from the Langanke effective potential; plain curves were computed from the KP II 
potential. Curve labels refer to reaction multipolarity and orbital angular momentum of the initial state. 




Ec.,„. (MeV) 

FIG. 8. Contributions of reaction mechanisms and incoming partial waves to the S factor, computed from the KP II potential. 
Cross sections were computed both using the plain Monte Carlo calculation, and normalizing the expected asymptotic form 
to the Monte Carlo operator densities; there are therefore two curves of each type, although they are indistinguishable for the 
total and for some of the mechanisms. Curve labels are as in Fig. lA 
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FIG. 9. The total El operator (solid curves) and the cross sections resulting from each of the smaller terms of the operator 
acting individually, shown as S factors. Plain curves were computed from the KP II potential. Curves with symbols x were 
computed from the Langanke potential. Labels for multipole and orbital angular momentum are as in Fig. M. Other labels 
are as in Eqs. (pl[-]46|), except that LWA2.1 and LWA2.2 denote the first and second terms of Eq. respectively. The data 
shown are the El cross sections and upper limits from Ref. [fl3|. 



and one from the KP II potential. The results are dominated by E2 transitions arising from D-wave scattering 
and by El transitions arising from P-wave scattering. Results from the two sets of scattering states are almost 
indistinguishable for these transitions, especially at low energies. In Fig. |^, the two sets of curves were both computed 
from the KP II potential. One set of curves was computed directly from the Monte Carlo calculation, while the other 
was computed from the expected asymptotic behavior of the integrand of Eq. [sj beyond 12 fm, as described above. 
Fig . ^ shows a further breakdown of the El operator into its various terms, along with the measured El contributions 
|l3| . It is seen that: 

1. ) The El cross section is larger than the data of Robertson et al. ||l^ by a factor of about 7 at 2 MeV. Most 
other models have also arrived at El cross sections that are too large relative to these data. A new aspect of our 
calculation is that the largest contribution is calculated from the nucleon-nucleon potentials as a relativistic center-of- 
mass correction. This corresponds closely to inserting measured cluster masses by hand in simpler models, as indicated 
by the fact that virtually identical (within 50%) results for the LWAc term are obtained from our six-body model 
and from a simple capture model using the same cluster potential, but with laboratory alpha and deuteron masses 
introduced. In Fig. ^, the top curve is the total El S factor (arising mainly from the center-of-energy correction), 
and the other curves indicate the S factors that would result from various small terms in the El operator individually. 
(These differ from their actual contributions to the S factor, since these contributions add coherently in the matrix 
element.) The center-of-energy and LWAb terms are obviously much more important than the other terms in our 
calculation, and the smaller terms do not affect the total significantly. 

We note also that Robertson et al. ||l^ report the opposite sign from what is expected for the El transition in their 
experiment. None of the effects in our calculation can produce an effect of the correct magnitude with the observed 
"wrong" sign, or reproduce the sharp drop in El amplitude from 1.63 to 1.33 MeV seen in the data. 

2. ) There is no significant Ml transition originating from the S'-wave scattering state at lowest order in LWA, and 
we estimate the next-to-lowest term to be much smaller than what we obtain for the lowest-order term. The Ml 
strength is expected to be generally very small, since it is zero by orthogonality arguments in a two-body model. 
However, at very low energy, any S'-wave contribution present may be expected to dominate the total cross section, 
and many-body wave functions can (in principle) provide small pieces of the final-state wave function which can be 
reached in such a transition. Such an effect is present in neutron captures on the deuteron and on '^He. The presence 
of a nonzero Ml contribution in our calculation raises two questions: to what extent is it a product of our choice 
of wave fimctions, and how will the presence of meson exchange currents (MEC), which is a factor-of-ten correction 
for the Ml cross section in '^IIe(n, 7y^e, affect our result? The answer to the first question lies in the discussion of 
orthogonalization at the end of Sec. ]l| and in Fig. orthogonal ground and scattering states produce much smaller 
Ml contributions than non-orthogonalized states, and the remaining Ml contribution is very sensitive to the choice 
of phcnomcnological scattering potentials, suggesting a near-cancellation that ought to be exact. The second question 
is less important than it may appear at first glance. The large relative MEC contributions to neutron captures 
on light nuclei typically arise from the isovector part of the Ml operator. The ad capture process is isoscalar; the 
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FIG. 10. Thermally- averaged reaction rates, for the Langanke (solid) and KP II (dashed) potentials. The upper panel shows 
rates scaled to the temperature dependence expected for a constant S factor (to make small differences visible), while the lower 
panel shows ratios to the rate of Robertson et al. jl^. Shown for comparison are the Robertson et al. rate (long-dashed), and 
the rate presently recommended in the NACRE compilation of reaction rates (dash-dotted) [^. The line at the top indicates 
the approximate temperature range for big-bang production of ®Li. 

isoscalar MEC operator is both significantly smaller than the corresponding isovector operator, and quantitatively less 
certain in magnitude. We have therefore neglected MEC contributions to the Ml operator. Only electric multipoles 
contribute to Coulomb breakup, so the discrepancy between our results and the indirect data (circles in Fig. |^) is not 
the result of omitting MEC. 

To compare our results with those obtained in a simpler model, we have also computed ad capture by treating 
the two clusters as point particles, and the ground state as an S'-wave bound state of these two clusters. We used 
the modified Langanke potential, described above, to compute both states. With this potential, the ^Li ground 
state energy is 0.07 MeV too high, and this affects the energy dependence of the cross section. The asymptotic 
normalization of the wave function, which sets the scale for the nonresonant cross section, depends in any model on 
the detailed short-range behavior of the ground state, and is larger for this model than for our six-body ^Li ground 
state. While the 711 keV resonance is equally well-reproduced in strength and width by both two-body and six-body 
models, the total cross sections away from resonance are as much as 16% smaller in the six-body model. The simple 
model also yields an El strength which is greater by about 30% at 3 MeV. The source of the greater El strength is an 
interplay between the differing asymptotic normalizations, the treatment of the center-of-energy correction (treated 
in the simple model by using laboratory values of alpha and deuteron masses) , and differences in the sizes of LWA2 
and LWAb terms. 

We have also calculated thermally-averaged reaction rates suitable for use in astrophysical calculations. We integrate 
from 0.1 keV to 10 MeV at temperatures from 10'' K to 3 x 10^° K to produce the quantity Na{(Jv) customarily used 
in reaction network calculations. Our result, shown in Fig. ^0|, is as much as 40% larger than previous estimates over 
a wide range of temperatures. The sources of these differences are not clear. They probably arise from several specific 
details of the capture models used, which have been normalized to the data using assumptions about the relative roles 
of the El and E2 components to produce the rates recommended in compilations. 

V. DISCUSSION 

The total cross sections we have computed compare favorably with the data in the region from the 3^ resonance to 
3 MeV. This reflects mainly a combination of reasonable fits of the scattering states to experimental phase shifts, and 
reasonable reproduction of the asymptotic behavior of the ^Li ground state in the S'-wave ad channel. Our calculation 
of the El cross section from the relativistic center-of-energy correction, and examination of several smaller corrections, 
is new; however, like all other models not actually scaled to the laboratory data, it is higher than the data. In our 
calculation, the El transition contributes just under 70% of the total cross section at 50 keV. 

A very large cross section for alpha-deuteron capture would be of some interest for cosmology, since it would allow 
significant production of ^Li in the big bang. However, we find that big-bang nucleosynthesis calculations using our 
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new cross sections can produce ^Li at a maximum level of 0.000 60 relative to ^Li, or ^Li/H=5 x lO"^**. The maximum 
occurs at a baryon/photon ratio of 2 x 10~^°, at the low end of the range of possible values suggested by other light- 
element abundances. At the higher baryon/photon ratio of 5 x 10~^" suggested by recent extragalactic deuterium 
abundance measurements, there is a factor of three less ^Li, but because of the increasing ^Li production, ^Li/^Li = 
0.000 07. Since the abundance ratio found in low-metallicity halo stars so far is in the vicinity of ^Li/'^Li = 0.05, it 
seems very unlikely that standard big-bang nucleosynthesis could account for observed abundances. It should be kept 
in mind that non-standard models not relying solely on the alpha-deuteron capture mechanism could, in principle, 
produce interesting amounts of ^Li ||3^ . 

In the future we expect to make several improvements on the present work. While the VMC wave functions for the 
s-shell nuclei give energies within 2% of the experimental value, the trial functions for p-shell nuclei are not as good. 
It should be worthwhile to try the more accurate GFMC ground states for the p-shell nuclei and calculate mixed 
estimates of the transition operators . Eventually one would also like to construct the scattering states with the 
GFMC method, but this will be considerably more difficult. Further, new three-body potentials that give a better fit 
to the energies of the light p-shell nuclei are currently under development and should improve both the VMC and 
GFMC results. An important aspect of these studies will continue to be the construction of ground states that have 
the proper asymptotic cluster behavior. This must be imposed in the variational trial function, because the GFMC 
algorithm is not sensitive to the tail of the wave function and will not be able to correct for any failures. 

The present calculation should serve as a useful starting point for studies of several other reactions of astrophysical 
interest involving p-shell nuclei, including ^He(i,7)^Li, ^He('^He, 7)^Be, and ^Be(p, 7)^6. Work on these reactions is 
now in progress. 
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